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Microquasars as sources of positron annihilation radiation 
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Abstract. We consider the production of positrons in microquasars, i.e. X-ray binary systems that exhibit jets 
frequently, but not continuously. We estimate the production rate of positrons in microquasars, both by simple 
energy considerations and in the framework of various proposed models. We then evaluate the collective emissivity 
of the annihilation radiation produced by Galactic microquasars and we find that it might constitute a substantial 
contribution to the annihilation flux measured by INTEGRAL/SPI. We also discuss the possible spatial distri- 
bution of Galactic microquasars, on the basis of the (scarce) available data and the resulting morphology of the 
flux received on Earth. Finally, we consider nearby "misaligned" microquasars, with jets occasionally hitting the 
atmosphere of the companion star; these would represent interesting point sources, for which we determine the 
annihilation flux and the corresponding light curve, as well as the line's spectral profile. We discuss the possibility 
of detection of such point sources by future instruments. 
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1. Introduction: Galactic Positron Sources 

Ever since the first detection of the gamma-ray signature 
(511 keV radiation) of the annihilating positrons in the 
Galaxy in 1970 (Johnson, Harnden and Haymes, 1972), 
the problem of the birth, propagation, interaction, and 
annihilation of the positrons has been a major topic of as- 
trophysical investigation, both observationally and theo- 
retically. Nonetheless, some of its fundamental issues, par- 
ticularly the questions of positron origin and propagation, 
have remained mostly unsolved puzzles. 

A series of observational campaigns, conducted by bal- 
loons 1 in the first two decades and then by satellites 2 , 
have produced increasingly more precise, high-resolution 
spectra. Recent data have allowed a determination of the 
physical conditions in the interstellar medium where the 
positrons annihilate (Guessoum et al. 2004; Churazov et 
al. 2005; Jean et al. 2006), based on an improved un- 
derstanding of the various physical processes involved 
(Guessoum, Jean & Gillard 2005). In parallel, consider- 
able progress has been made in mapping the Galaxy at 
511 keV, first by CGRO-OSSE (Purcell et al. 1994; Cheng 

the pioneering Rice University group's campaign, the Bell- 
Sandia group's, GRIS, HEXAGONE 
2 HEAO 3, CGRO-OSSE, WIND-TGRS, INTEGRAL-SPI 



et al. 1997; Purcell et al. 1997; Milne et al. 2000; Milne et 
al. 2001) and, in particular, recently by INTEGRAL-SPI 
(Knodlseder et al. 2005). 

Despite that substantial progress, the origin of the 
huge amounts of positrons produced in the Galaxy (10 43 
e + s _1 ) is still eluding us. A large number of potential 
sources has been proposed over the years: cosmic ray in- 
teractions with the interstellar medium (Ramaty, Stecker 
& Misra 1970), pulsars (Sturrock 1971), radioactive nu- 
clei produced by explosive nucleosynthesis in supernovae 
(Clayton 1973) or novae (Clayton & Hoyle 1974), com- 
pact objects housing either neutron stars or black holes 
(Ramaty & Lingenfelter 1979), matter expelled by red gi- 
ants (Norgaard 1980) and Wolf-Rayet stars (Dearborn & 
Blake 1985, Prantzos and Casse 1986), gamma-ray bursts 
(Lingenfelter & Hueter 1984), (light) dark matter (Rudaz 
& Stecker 1988; Boehm et al. 2004), low-mass X-ray bina- 
ries (Prantzos 2004), hypernovae (Casse et al. 2004), mil- 
lisecond pulsars (Wang, Pun, & Cheng 2006), and, most 
recently, pair production from the collision of gamma-ray 
photons from "SMall Mass Black Holes" (SMMBHs) and 
X-ray photons from the Galactic Center black hole Sgr A* 
(Titarchuk & Chardonnet 2006). 

The recent mapping of the Galaxy at 511 keV 
(Knodlseder et al. 2005), particularly the very large 
Bulge/Disc (B/D) ratio inferred, has placed rather se- 
vere constraints on the positron candidate sources. Of 
the aforementioned potential sources, Type la supernovae, 
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low-mass X-ray binaries (LMXB's), and dark matter have 
emerged, for various reasons, as possible (or, rather, not 
impossible) candidates. However, each of these has its own 
difficulties: 

1) Although the total positron emissivity of SNIa in 
the Galaxy matches the observed one, the corresponding 
bulge emissivity is estimated (e.g. Prantzos 2004; Diehl, 
Prantzos & von Ballmoos 2005) to be an order of mag- 
nitude lower than required by SPI observations. Even as- 
suming that large systematic uncertainties affect that esti- 
mate, the expected galactic distribution of SNIa docs not 
correspond to the SPI B/D ratio, unless additional as- 
sumptions are made (e.g. important number of currently 
undetected SNIa in the galactic bulge, or transfer of a 
large fraction of the disk positrons to the bulge through 
the galactic magnetic field, see Prantzos 2006); 

2) Low-mass X-ray binaries (LMXBs) display a spatial 
distribution considerably peaked toward the Galactic cen- 
tre region (Grimm, Gilfanov & Sunyaev 2002), and their 
global energetics allow for a substantial positron produc- 
tion galaxywide (Prantzos 2004). However, the brightest 
of them (in X-rays) lie in the disc, not in the bulge of 
the Milky Way (Prantzos 2004); if positron emissivity cor- 
relates with X-ray emissivity (as is commonly assumed), 
LMXB's should be excluded on the basis of inadequate 
morphology. 

3) Dark matter particles as sources of Galactic 
positrons may also suffer from "morphology" problems, or 
at least uncertainties, but also from inconsistencies in the 
framework of particle physics (see Ascasibar et al. 2005). 

Finally, we must note the important constraint that 
Beacom & Yiikscl (2005) have pointed out w.r.t. all mod- 
els in principle: if the positrons are relativistic (E ^ 3 
MeV) , then the annihilation in flight of a fraction of them 
(which is very large when the medium is neutral) will pro- 
duce continuum photons of high energies (> a few MeV) 
that could be detected by instruments such INTEGRAL 
and COMPTEL. This result places strong constraints on 
models, particularly those that have positrons produced at 
high energies and annihilating "in flight" , i.e. while slow- 
ing down. 

In the case of LMXBs, positrons should be produced 
as e + -c~ pairs in the inner regions of their accretion 
discs. Some of these would annihilate locally, but a non- 
negligible fraction would be channeled out by jets — when 
these exist. In general, positrons would reach the ISM rela- 
tively far from their source, and they would propagate and 
annihilate, contributing to the diffuse 511 keV emission. 
In special cases, when the jets are strongly inclined ("mis- 
aligned" ) toward the plane of the binary system, positrons 
could periodically hit the atmosphere of the companion 
star where they would produce an annihilation signature, 
characterized, in particular, by its time variability and line 
profile. 

In this paper we consider the possible fate of positrons 
ejected by microquasar jets. We first estimate the rate 
of positron production in jets, by reviewing estimates 
from other studies and by using simple energetics con- 



siderations; in that respect, we pay special attention to 
upper limits derived from the INTEGRAL-SPI measure- 
ments of positron annihilation fluxes from point sources in 
the Galaxy. We then consider two possible consequences: 

1) the collective microquasar contribution to the flux of 
galactic annihilation radiation, which we find to be poten- 
tially substantial (especially from the central regions); and 

2) the annihilation line signature of "misaligned" micro- 
quasars. In the latter case, we study the fate of positrons 
in the companion's atmosphere in order to compute the 
photon emissivity, and determine the characteristic "light 
curve" of such a source; we find that interesting obser- 
vational signatures could be obtained from close enough 
and/or strongly active sources. 

In Sect. 2 we present a brief overview of X-ray bina- 
ries, microquasars, and jets. In Sect. 3 we present our es- 
timates of the rate of positrons produced by microquasar 
jets. In Sect. 4 we consider the overall contribution of mi- 
croquasars to the global galactic annihilation flux. In Sect. 

5 we consider the annihilation of positrons in misaligned 
microquasars. In Sect. 6 we summarize our conclusions. 

2. X-Ray Binaries, Microquasars, and Jets 

X-ray binaries (XRB's) are systems containing a com- 
pact object (cither a neutron star or a stellar-mass black 
hole) accreting matter from a companion star. A total 
of 280 galactic X-ray binaries are currently known: 131 
High Mass X-ray Binaries (HMXB's) (Liu, van Paradijs 

6 van den Heuvel 2000), in which the companion star has 
a mass > 5 solar masses (spectral type O or B) and where 
the mass transfer usually takes place by way of the strong 
stellar wind; and 149 Low Mass X-ray Binaries (LMXB's) 
(Liu, van Paradijs & van den Heuvel 2001), in which the 
companion has a low mass (spectral type later than B), 
and the mass transfer is carried out by Roche lobe over- 
flow. Grimm, Gilfanov & Sunyaev (2002) have shown that 
HMXB's tend to be distributed along the galactic plane, 
while LMXB's tend to be clustered in low Galactic longi- 
tudes. The number of XRB's brighter than 2 x 10 34 erg/s 
is estimated to ^700 in the Galaxy (Grimm, Gilfanov & 
Sunyaev 2002; Paredes 2005). 

Among the detected XRB's, 43 have now been found to 
exhibit radio emission, which is interpreted as synchrotron 
radiation. It is generally believed that the radio emission 
is evidence of jets. Of the 43 objects, 35 are LMXB's and 
8 are HMXB's. Moreover, of the 43 systems, 16 are con- 
firmed cases with resolved jets (Ribo 2005, Paredes 2005, 
Chaty 2005). Such systems (a compact object accreting 
from a companion star and ejecting a stream of relativis- 
tic particles) are known as microquasars. It is likely that 
all 43 radio-emitting X-ray binaries (REXB) are micro- 
quasars. The population of microquasars in the Galaxy is 
estimated at about 100 (Paredes 2005), although only 16 
are "officially" known (or confirmed) at present. We note 
that Pandey et al. (2006a, 2006b) have recently reported 
4 new REX sources: IGR J17091-3624, IGR J17303-0601, 
IGR J17464-3213, and IGR J18406-0539; in addition, two 



N. Guessoum et al.: Microquasars as sources of positron annihilation radiation 



3 



sources, X Nor X-l and IGR J17418-1212, have been pre- 
sented as microquasars (Klein- Wolt, Homan & van der 
Klis 2004, and Tsarevsky 2004, respectively). We refer to 
these six sources as "new tentative microquasars". The 
catalog of microquasars is expected to grow rapidly with 
the array of instruments and studies now investigating 
these objects (from radio to gamma rays). 

Although jets of such systems were first observed in 
1979 (in the peculiar object SS 433), microquasars were 
identified and imaged in the Galaxy only in 1992 when 
Mirabel et al. (1992) performed high-resolution radio ob- 
servations of the GC's "great annihilator" IE 1740.7 - 
2942 and soon afterwards in GRS 1915+105 (Mirabel & 
Rodriguez, 1994). 

Recent studies (Falcke, Kording & Markoff 2004; 
Mirabel 2004; Fender, Belloni & Gallo 2004, 2005) consid- 
erably improved our understanding of the conditions that 
lead to the emergence of jets in XRBs. First and fore- 
most, a connection has been established between X-ray 
flux and jet formation: jets seem to appear when the ac- 
cretion disc X-ray luminosity is low. Jets are apparently 
produced when the inner disc is replenished; there is a 
clear general pattern of steady jets in the "low/hard state" 
of the X-ray (microquasar) sources, while no jet is seen in 
"the high/soft state". 

A relatively elaborate model of this scenario has been 
proposed by Fender, Belloni & Gallo (2004, 2005). A 
relation between the power in the jet (when it exists) 
and the X-ray luminosity of a given source is derived: 
Lj et = A steady L^x with luminosities expressed in units 
of the corresponding Eddington luminosity. These authors 
argue that steady jets are produced when the X-ray spec- 
trum of a source hardens beyond a certain value (which 
may be universal or vary somewhat from one source to 
another). The spectrum softens when the X-ray luminos- 
ity increases above about 1 % of the Eddington value; 
when this happens, the jet first increases in speed and 
then quickly gets suppressed and disappears. 

There seems to be an agreement on the above rela- 
tion between L je t and L^ 5 , which is supposed to hold for 
sources individually. On the other hand, there seems to be 
considerable uncertainty over the value of A stea dy (Fender, 
Maccarone & van Kesteren, 2005); indeed, using various 
sources (e.g. XTE J1118+480) for "calibration", authors 
obtain values ranging from 0.006 as a lower limit (Fender, 
Belloni & Gallo 2005) to 0.3 (Malzac, Merloni & Fabian, 
2004), which Fender, Belloni & Gallo (2005) take as an up- 
per limit. In the case of transient jets, Fender, Maccarone 
& van Kesteren (2005) conclude that the value of A tra ns 
would range between 0.04 and 4.0, at least for black hole 
sources, which tend to exhibit jets that are about 10 times 
more powerful than those of neutron star jets, both in the 
steady and in the transient cases. 

Secondly, a general correlation is also found between 
the velocity of the outflow and the X-ray luminosity of 
the accreting source: increases in X-ray luminosity tend 
to accelerate the jets as long as the source remains in the 
low/hard state. Fender, Belloni & Gallo (2005) further 



argue that the velocities of transient jets are significantly 
larger ( Z, 0.87 c) than those of steady jets ( ;$ 0.7 c). Jets 
with larger Lorentz factors have also been considered, but 
the standard cases are those presented by Fender, Belloni 
& Gallo (2005). 

The particle content of the jets is among the major 
unresolved issues in the field of microquasars (e.g. Schcck 
et al. 2002). Observations of emission lines in the case 
of SS 433 (Marshall, Canizares & Schulz 2002) suggest 
a substantial baryonic content of its jets. Several studies 
have considered the implications of a hadronic jet compo- 
sition, either for the galactic cosmic ray content (Heinz & 
Sunyaev 2002, Fender, Maccarone & van Kesteren 2005) 
or for the synthesis of Li on the surface of the compan- 
ion star (Butt, Maccarone & Prantzos 2003). Others have 
considered the gamma-ray and neutrino production at the 
surface of the companion by impinging high-energy pro- 
tons (Romero & Orellana 2005). On the other hand, ar- 
guments for e + -e~ pair-dominated jets have been put for- 
ward, especially in cases involving the extraction of the 
spin energy of the black hole (Celotti & Blandford 2001). 
A strong argument in favor of such a leptonic composi- 
tion is the repeated observations of highly polarized jets. 
Several leptonic microquasar models have been proposed 
in recent times, e.g. Bosch-Ramon and co-workers (2004a, 
2004b, 2005), Dermer & Bottcher (2006); see the review 
by Romero (2005). Some rather complex models have also 
been proposed, like "two-flow" models (pair beam sur- 
rounded by a mildly relativistic e~-p plasma) , or e~-p 
jets that later get loaded with pairs by interactions with 
high-energy photons (e.g. Scheck et al. 2002 and references 
therein). 

In this work we consider microquasar jets channel- 
ing positrons (or c + -c~ pairs) from the inner regions of 
the source's accretion disc into the ISM or towards the 
companion star (if the jet is "misaligned"). We also note 
that positrons can also be produced by hadrons from the 
jet colliding with nuclei (e.g. through p-p^p + n + 
7r + ), but such processes would contribute only negligible 
amounts of positrons due to the small numbers of high- 
energy hadrons as well as the low values of the relevant 
cross sections. 

In order to estimate the flux of the resulting annihi- 
lation radiation, the rate of positron ejection by the jets 
must then be determined. 

3. Rate of positrons produced by microquasar 
jets 

The high-temperature, high-density conditions in the in- 
ner regions of the binary system's accretion disc produce 
electron-positron pairs via 7 + 7 — > e + + e~ reactions. 
A fraction of them annihilate close to the compact object, 
but when jets appear, they may channel out a significant 
number of pairs. A few authors have attempted to model 
and estimate the production and ejection of pairs in those 
conditions. 
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According to Beloborodov (1999) the pairs are cooled 
to energies of 1 - 10 keV and blown away by soft radi- 
ation to form a semi-relativistic wind. Depending on the 
compactness of the source (I = Lot / 'm e (? R, where L and 
R are the power and radius of the emitting region, and 
&t is the Thomson cross section), the plasma will form 
an optically thin or thick atmosphere. The density of the 
outflow and the rate of ejection of electrons and positrons 
depend on the rate of production of pairs (which, in its 
turn, depend on the photon "seed" spectrum and on the 
accretion disc model), the annihilation rate and the "es- 
cape efficiency" of the pairs. Under the assumption of an 
optically thin pair wind, where the pairs escape before 
they annihilate {r esc < r anni with r ann ~ l/n e u T c, where 
n e is the electron density), Beloborodov (1999) shows that 
the maximum pair luminosity is given by 



which translates into a rate of pair injection in the jet of 
~4x 10 41 s _1 . A substantial fraction of these pairs, per- 
haps up to 90 % (as argued by Misra & Melia 1993), would 
annihilate near the base of the jet, producing a broad and 
rcdshifted line (that is rather difficult to detect); as many 
as ~ 10 41 s _1 positrons are then expelled into the ISM 
generally or, occasionally, in the direction of the compan- 
ion star. 

Misra & Melia (1993) suggest that the intense radia- 
tion field is responsible for the Compton acceleration of 
the pairs produced in the inner regions of the accretion 
disc. They find that a large rate of pairs (up to 6 x 10 42 
s _1 ) stream outwards from the disc (at velocities of ~ 0.7 
c), even after 90 % have annihilated near the base. One 
must note, however, that this large rate is obtained with 
an accretion rate of ~ 5 x 1O~ 8 -M yr _1 , a rate that can 
normally be attained only in episodic outbursts. Indeed, 
the model of Misra & Melia (1993) was mainly aiming 
to reproduce the IE 1740.7 -2942 "annihilation flare" of 
1991. 

Yamasaki, Takahara & Kusunose (1999) consider two- 
temperature accretion discs by taking into account the 
formation of rclativistic pair outflows, in both AGN and 
microquasars. They show that in the inner regions of the 
discs, when the mass accretion rate becomes larger than 
about one tenth of the Eddington rate (1.4 x 1O 17 M/M 
g s _1 where M is the mass of the compact object) or ~ 
2 x 1O~ 1O M yr _1 , most of the viscously dissipated energy 
is converted into the thermal and kinetic energy of the 
electron-positron pairs. They obtain a maximum power 
of the pair outflow of 0.136 Lsdd for & n accretion rate 
of 1O~ 9 M yr- 1 (assuming M = 1M and R = 10i? o ), 
which translates into a pair ejection rate of 2 x 10 42 s~ 4 , 
assuming the jet is leptonic. 

Considering the global energetics of microquasar jets 
in the Galaxy, one may note that 

- at L x = 0.5 L Edd one obtains N e + = 2 x 10 42 s _1 
(Yamasaki, Takahara & Kusunose 1999), which in the 
following estimate we may use as a "yardstick" ; 



— steady jets are produced at Lx = 0.01 - 0.1 L^dd 
(Fender, Bclloni & Gallo 2004, 2005); 

and 

— N e + cx Lj et oc L x , where a=0.5 or 1 (see discussion 
in Sec. 2) and where the value of the proportionality 
constant A is not needed if one uses the above "yard- 
stick" ; 

one can infer that the positron emissivity of a jet in the 
steady state lies in the range N e + ~ 4 x 10 40 — 9 x 10 41 
c+ s _1 , with 10 41 e + s _1 as a reasonable average value. 

Finally, in attempting to fit the Nova Muscae 1991 "an- 
nihilation flare" , Kaiser & Hannikainen (2002) inferred a 
positron annihilation rate of 2 x 10 43 s _1 in the atmosphere 
of the companion star (assuming a distance to Earth of 
w 5.5 kpc). They argued that such a high pair yield 
can be achieved when the plasma is "photon-starved", 
that is when the number of high-energy photons in the 
disc strongly exceeds that of the soft seed photons (see 
Zdziarski, Coppi & Lamb 1990). 

Despite the wide range in the results obtained in the 
aforementioned studies, one concludes that a reasonable 
estimate for "steady-state" production and ejection of 
pairs in normal conditions is ~ 10 41 s _1 . 

Taking a different, empirical approach, interesting 
upper-limit constraints to the microquasar positron emis- 
sivity may be established by way of the recent annihilation 
radiation measurements by INTEGRAL-SPI. Knodlcseder 
et al. (2005) have published 3-cr flux upper limits for a 
dozen galactic microquasars/LMXB's/HMXB's; we have 
completed the data for other sources of interest to us here, 
particularly the "misaligned" microquasars that we will 
consider in some detail in Sec. 5. Table 1 lists the sources 
we have considered, with the corresponding SPI upper- 
limits on their 511 keV emission flux (taking the sources 
to be point-like) and inferred rates of positron injection, 
assuming that positrons do not annihilate far from the 
source. The bottom part of the table lists microquasar 
sources, the last two or three being misaligned ones. 

It can be seen that the data of the Table imply that 
upper limits for the steady positron annihilation rate of in- 
dividual sources are always >10 41 e+s" 1 (albeit with large 
uncertainties, due in particular to uncertain distance esti- 
mates) , so the upper-limit rates inferred from INTEGRAL 
data do not clash with the average value derived in the 
works cited above. 

From the different considerations and estimates pre- 
sented above, it appears that a "canonical average" rate 
of ~ 10 41 c + s _1 ejected by the microquasar jets is not un- 
reasonable. We will adopt this value for the determination 
of expected fluxes from individual sources (section 5), but 
realizing the corresponding large uncertainty, we will also 
consider lower (10 40 e + s _1 ) values for "weak-episode jets" 
as well as larger (10 42 c+s -1 ) ones, in "strong flare" cases. 
We immediately note that should this average, "canonical 
value" prove to be reasonably correct, the ^100 micro- 
quasars that are believed to exist in the Galaxy (Paredes 
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Table 1. Limits on positron rates from SPI upper limits for XRB sources of interest. The positron rates were calculated 
assuming a positronium fraction of 95%. The bottom part of the table lists microquasar sources, the last two or three 
being misaligned ones. 



Source 


Type 


Distance 


1 


b 


3<r Flux Limit 


Positron Rate 






(kpc) 


(deg) 


(deg) 


(10~ 4 cmV 




(e+s- 1 ) 


GX 349+2 


LMXB 


Z 10 


349.1 


2.75 


U.o 




"••> -| mm I ,-,49 

> 1.7 x 10 


GX 5-1 


LMXB 


8 


5.08 


-1.02 


U. 1 




9.4 x 10 41 


Nova Muscae 


LMXB 


3 


295.3 


-7.07 


o 
1 




3.8 x 10 41 


A 0620-00 


LMXB 


2 


209.96 


-6.54 


O.O 




3.1 x 10 41 


Cen X-4 


LMXB 


1.2 


332.24 


23.89 


i 

1. 1 




5.2 x 10 


GRS 1915+105 


LMXB 


12.5 


45.37 


-0.22 


1 




3.3 x 10 42 


Cir X-l 


LMXB 


10 


322.12 


0.04 


1.1 




2.3 x 10 42 


Cyg X-3 


HMXB 


9 


79.85 


0.7 


1 




1.7 x 10 42 


IE 1740.7-2942 


LMXB 


8.5 


359.1 


-0.11 


0.9 




1.4 x 10 42 


GRS 1758-258 


LMXB 


8.5 


4.51 


-1.36 


0.7 




1.1 x 10 42 


GX 339 


LMXB 


^8 


0.68 


-0.22 


0.8 




< 10 42 


SS 433 


HMXB 


4.8 


39.69 


-2.24 


0.9 




4.5 x 10 41 


LS 5039 


HMXB 


2.9 


16.88 


-1.29 


0.9 




1.5 x 10 41 


Sco X-l 


LMXB 


2.8 


359.1 


23.78 


1.5 




2.4 x 10 41 


Cyg X-l 


HMXB 


2.5 


71.33 


3.07 


1 




1.4 x 10 41 


XTE J1118+480 


LMXB 


2.5 


157.66 


62.32 


4.5 




5.9 x 10 41 


LS I +61° 303 


HMXB 


2 


135.68 


1.09 


3.3 




2.8 x 10 41 


IGR J17091-3624 


? 


8.5? 


-10.48 


2.21 


0.7 




1.1 x 10 42 


IGR J17303-0601 


LMXB 


8.5? 


17.93 


-1.61 


0.9 




1.4 x 10 42 


IGR J17464-3213 


LMXB 


8.5? 


-2.87 


15.01 


0.9 




1.4 x 10 42 


IGR J18406-0539 


? 


8.5? 


26.67 


-0.17 


1 




1.5 x 10 42 


X Nor X-l 


LMXB 


8.5? 


-23.09 


0.25 


1 




1.6 x 10 42 


IGR J17418-1212 


? 


8.5? 


13.93 


9.41 


0.9 




1.4 x 10 42 


XLE J1550-564 


Misaligned? 


5.3 


-34.12 


-1.83 


1.1 




6.6 x 10 41 


V4641 Sgr 


Misaligned 


9.6 


6.77 


-4.79 


0.7 




1.4 x 10 42 


GRO J1655-40 


Misaligned 


3.2 


-15.02 


2.46 


0.8 




1.7 x 10 41 



2005) would produce a global annihilation cmissivity near 
that measured by SPI (and previous instruments). 

Another upper limit on the collective emission by 
galactic microquasars may be derived by considering 
the global energetics. Indeed, (i) the total luminosity 
of LMXRBs in the Milky Way is 2 - 3 10 39 erg s" 1 
(Grimm, Gilfanov & Sunyaev 2002), while (ii) the lumi- 
nosity of - 10 43 e+ s" 1 observed by SPI/INTEGRAL is 
~ 10 erg s _1 , assuming the positronic jets arc mildly 
relativistic (i.e. positrons have energies <1 McV). Of 
course, as Grimm, Gilfanov & Sunyaev (2002) note, the 
Galactic LMXRB population is dominated by a dozen 
bright sources (of ~10 38 erg s _1 each) lying in the disk. 
But the remaining fraction of ~ 90 % is highly clus- 
tered towards the bulge (Fig. 1 in Grimm, Gilfanov & 
Sunyaev 2002), as required by the SPI data and noted in 
Prantzos (2004), while their collective luminosity is 2 - 3 
10 38 erg s _1 , i.e. 20 times larger than required to explain 
the Galactic positron energetics. We note that those are 
precisely the low-luminosity (sub-Eddington) sources that 
may produce jets (see the discussion in Sec. 2). According 
to estimates by Paredes (2005), based on current micro- 
quasar statistics, there are about 100 microquasars in 
the Milky Way, which corresponds to about 1/3 of the 



~ 300 LMXRBs in our Galaxy as estimated by Grimm, 
Gilfanov & Sunyaev (2002). Applying this correction fac- 
tor of 1/3 to the microquasar energetics, still leaves about 
10 38 erg s _1 available for positron production in their jets, 
i.e. about 6 times more than required by SPI data. We 
note at this point that some of the currently observed 
microquasars result from HMXRBs, but their fraction is 
rather small (less than 20 %) and does not affect the ar- 
guments presented here (see next section for a more de- 
tailed treatment). Moreover, we wish to stress the large 
(order-of-magnitude) uncertainty that exists in our cur- 
rent knowledge of the ratio between the energy that goes 
into positrons or into the jet (as a whole) compared to the 
energy t hat is radiated, an issue that is further compli- 
cated by the unknown lepton-to-hadron content ratio of 
the jet. 

The simple estimate made in the previous paragraph 
implies that the ratio between the positron power and the 
X-ray luminosity of microquasars is, on average, about 
16 %. One must recall, however, that there is a relation 

1/2 

between the two quantities, namely N e + ~ Lj et ~ L^ 
(references given above), which leads to more realistic ex- 
pectations. 
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Indeed, calculating the derivative of the cumulative 
X-ray luminosity function given by Grimm, Gilfanov & 
Sunyaev (2002)' for LMXRB's (see Eq. 15 of Grimm, 
Gilfanov & Sunyaev 2002) and assuming that 1/3 of 
LMXRBs arc microquasar, we derived a differential lu- 
minosity distribution function for microquasars: 



dN/dL 



x 



2.6L 



c 1 ' 26 , (2) 

where Lx is in Eddington units (1.3xl0 38 ergs/s) and 
relates to L e + as L e + = B L x , with a = 1 for a linear 
relationship and a = 1/2 for the alternate, more realistic 
relation. 

Now if the total positron power in the Galaxy is re- 
quired to be 10 37 erg/s (i.e. 0.1 in Eddington units), and 
using the above equation (Eq. 2) to integrate for the total 
power, one finds that for a = 0.5, B is between 0.9 and 
1.6 % depending on the assumption made on the maxi- 
mum X-ray luminosity of the microquasar (1 or 0.1 I^Edd, 
respectively), which are very reasonable figures. 

Furthermore, having estimated B and thus obtained 
an analytic distribution function for the microquasars' 
positron emissivity, one can obtain the total rate of 
positrons emitted by all microquasars (of different lumi- 
nosities) in the Galaxy by integrating dN/dLx against 
the luminosities and the rate of positrons emitted by each 
source, as obtained by Yamasaki, Takahara & Kusunose 
(1999). Again, for the linear relation (a = 1) one obtains a 
total rate N e + itot between 0.3 and 1.4 x 10 43 e + /s (depend- 
ing on the maximum X-ray luminosity of a given micro- 
quasar, i.e. 0.1 or 1 I^Edd, respectively); however, for the 
preferred non-linear relationship (a = 1/2), one obtains 
N e + itot between 1.8 and 3.1 x 10 43 e + /s (microquasar X- 
ray luminosities of 0.1 and 1 l^Edd, respectively), values 
that are very encouraging in considering microquasars as 
possible contributors to the overall positron annihilation 
flux from the central galactic regions. 

4. Contribution of microquasars to the galactic 
positron annhilation radiation 

Assuming that ~ 100 microquasars exist in the Galaxy 
(Paredes 2005), and noting from the locations plot of the 
22 microquasars presently known/confirmed (Figure 1) 
that about half of them are in the central regions, i.e. 
~ ±25 degrees from the GC, we can estimate the flux 
of annihilation radiation that can be expected from such 
a population of sources with jets ejecting w 10 41 e + s _1 
on average and compare it to the global galactic centre 
annihilation flux reported by Knodlseder et al. (2005). 

The microquasars can have either "misaligned" (highly 
inclined with respect to the orbital plane) or "normal" 
(low inclination) jets. The normal microquasars will pour 
out their positrons into the ISM, where the usual 2-3 pho- 
ton production processes take place (see Guessoum, Jean, 
& Gillard 2005). In the case of misaligned microquasars 
(about 15-20 % of the total) we assume, following Butt, 
Maccarone & Prantzos (2003), that their jets hit the com- 
panion ~10 % of the time. Only 1 of the ~ 2 annihilation 



photons will emerge then, the other being absorbed in the 
companion's atmosphere (see Sec. 5). The remaining 90% 
of the positrons are thrown into the ISM, where they pro- 
duce 2-3 photons each. 

The total flux of annihilating positrons coming from 
the inner Galaxy would be the sum of the flux contributed 
by the misaligned jets and the normal jets: 



F 



where 



and 



Fyriia 



+ F n . 



4ttD 2 



X 2 X /line 



(3) 



(4) 



f/j.Q A M Q ,/mis^Q N e 

AttD 2 

x[/jet-* + 2 X (1 - 



/jet — *)] X /line , 



(5) 



where N^q is the total number of microquasars believed to 
exist in the Galaxy (~ 100), J^q is the fraction of ^Q's we 
assume to be in the inner regions of the Galaxy (« 50%), 
/mis/xQ is the fraction of /zQ's assumed to have misaligned 
jets (« 20%), /jot-* is the fraction of the time a misaligned 
jet will hit the companion star's atmosphere (« 10%), N e + 
is the rate of ejection of positrons from a typical jet (10 41 
s _1 ), /ii no is the fraction of photons emitted with the line 
energy (511 keV) as opposed to continuum (0 - 511 keV) 
energies, and D is the distance to the Galactic centre (8.5 
kpc); /ii n o is obtained from ^/p s + (1 — /p s ), where /p s is 
the "Positronium fraction" , i.e. the fraction of positrons 
that annihilate via formation of Positronium (the bound 
e + -e _ system), which has repeatedly been found in galac- 
tic annihilation radiation measurements to be « 0.95 (see 
references given in Sec. 1); 

With those parameter values, we obtain: 



4 x 10~ 4 ph cm" 2 s" 



6 x 10~ 5 ph cm 



(6) 
(7) 



The total flux F w 5 x 10~ 4 ph cm~ 2 s _1 is within 
a factor of 2 of the SPI-measured annihilation flux (1.1 x 
10 ph cm _2 s _1 ), a result that is quite encouraging, con- 
sidering the uncertainties on different parts of the problem 
(mostly due to our currently limited knowledge of micro- 
quasar jet energetics). 

As noted in Sec. 1, the spatial distribution of the 
511 keV flux detected by SPI-INTEGRAL puts severe 
constraints on candidate sources of positrons. Although 
far from complete at present, the available sample of 
known/confirmed microquasars appears encouraging in 
that respect. Figure 1 displays the position and type 
(LMXB/HMXB) of the 22 currently known microquasars. 

While fully aware that this list represents only about 
one fifth of the microquasars believed to exist in the 
Galaxy, we can still attempt to determine a Bulge-to-Disc 
(B/D) ratio of the annihilation produced by such sources 
and compare that with the limits obtained from SPI data 
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(Knodlseder et al. 2005), which inferred a rate of positron 
annihilation of (1.5 ± 0.1) x 10 43 e + s~ 4 in the bulge and 
(0.3 ± 0.2) x 10 43 e+s" 1 in the disc. 

Referring to Table 1, which lists the types, positions 
and distances of these microquasars, we note that roughly 
80 % of them are LMXB's, about 15 % of which are in the 
halo, 35 % are in the disc, and about 50 % are in the bulge, 
while 20 % of the sources are HMXB's, ~ 10% of which 
may be in the bulge and the rest in the disc. The canonical 
average rate of positron production has been taken to be 
about 10 41 e + s _1 by a typical jet from an LMXB micro- 
quasar, while for HMXB microquasars jets are ten times 
less powerful and thus about three times less productive 
in positrons (recall that L e + e - ~ L- }ct ~ L^ 5 ). One must 
then take into consideration the confinement probability 
of positrons ejected from these sources; we note that: i) 
according to Jean et al. (2006), positrons produced in the 
bulge do not escape, they end up annihilating in the bulge 
if their energy is below ^10 MeV, which is the case for 
positrons produced by microquasars, if one ignores (as in 
this basic treatment) internal acceleration processes; ii) 
the scale height of LMXBs in the disk is ~400 pc, while 
the gas has a scale height of ^100 pc, so about 50 % 
of positrons produced by LMXB's there are ejected to- 
ward the disk and annihilate, while the rest (50 %) of the 
positrons are released in the halo and either contribute to 
a diffuse annihilation emission (unseen by spectrometers) 
or propagate following the galactic magnetic field lines to- 
ward the bulge (Prantzos 2006); hi) the escape fraction 
of positrons produced by HMXB's in the disk is uncon- 
strained, and for simplicity we take it to be the same as 
that of positrons from disk LMXB's. With these fractions, 
the net rate of positrons annihilating in the bulge is found 
to be « 4.1 x 10 42 e + s~ 4 (about one third the SPI bulge 
rate), while the net rate of positrons annihilating in the 
disc is found to be w 1.7 x 10 42 e + s~ 4 ; this would give 
a B/D ratio of 2.4, which is somewhat smaller than the 
lower SPI limit (B/D) spi > 3, indicating that positrons 
produced in the disk are escaping in large(r) fractions than 
assumed here. 

One can also turn these estimates around and set limits 
on positron production rates from microquasars, knowing 
that the net rates in the bulge and in the disc cannot ex- 
ceed those set by SPI as mentioned above. (Note that since 
we are taking a canonical average value of positron pro- 
duction rate for all LMXB microquasars and assuming the 
HMXB microquasars to produce only one third as many 
positrons, for luminosity reasons as explained above, the 
B/D ratio is independent of the jet positron rate and can- 
not help set constraints on it, at least in our simple model.) 
Assuming the ratios of LMXB /HMXB microquasars and 
the escape fractions given in the previous paragraph, we 
find that the SPI limits would be violated if positrons are 
produced at steady rates greater than ~ 3 x 10 41 c+s -1 . 
Furthermore, we note that our model easily satisfies the 
Beacom & Yiiksel (2005) constraint since our positrons 
have kinetic energies between 2.5 keV and 660 keV (jet 
speeds range from 0.1 c to 0.9 c), so their annihilation in 
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Fig. 1. Microquasars (position and type) in the Galaxy. 

flight does not produce continuum emission of photons at 
higher energies. 

Finally, we note that the line profile that can be ex- 
pected in this scenario is almost exactly the same as the 
one obtained by SPI and analysed recently (Churazov et 
al. 2005, Jean et al. 2006), since the positrons are here 
poured mostly into the ISM; the small fraction ejected 
from misaligned microquasars that annihilates on the 
companion star's surface will produce a line of similar 
profile, as the physical conditions (temperature and ion- 
ization) in the atmospheres of those stars are essentially 
identical (T ~ 8000 K, H/He composition, partially ion- 
ized) as those in the ISM phases where the positrons an- 
nihilate on large scales. 

5. Annihilation in the Atmosphere of the 
Companion Star 

In this section we consider in some detail the fate of the jet 
positrons hitting the atmosphere of the companion star, in 
the case of misaligned microquasars. Although this hap- 
pens in a small number of cases (~15-20% of all the micro- 
quasars, see Sec. 4) and for a small fraction of the XRB's 
period (~10%) there are potentially interesting gamma 
ray signatures of such an interaction. 

Positrons arrive in the companion star with speeds 
ranging from 0.1 c to 0.9 c (canonical value Vj et ~ c/2), i.e. 
kinetic energies of 2.5 keV to 660 keV. They lose energy 
in the stellar atmosphere which has a very sharp density 
profile (see Fig. 1 in Guessoum & Jean 2002). Using stellar 
atmosphere models presented in Guessoum & Jean (2002) 
we estimate the distance for the energetic positrons to 
slow down below 60 eV to be A_e ^ 0.1 g cm -2 . At this 
energy, positrons are able to form a positronium by charge 
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exchange with H atoms. Low energy positrons that do not 
form a positronium quickly thermalize and annihilate in- 
situ. The maximum atmospheric depth at which positrons 
annihilate (Xe) is lower than the mean free path of 511 
keV photons (A 511kcV ~6.6 g cm -2 ) in the atmosphere. 
Then in these conditions, most 0.95%) of the anni- 
hilation photons that are emitted upward do escape the 
atmosphere of the companion star. 

This assumption can be considered valid as long as 
the jet impact does not strongly deform the stellar at- 
mosphere. It should be noted, however, that the energy 
deposited by the jets into the star's atmosphere (~ 10 36 
erg s -1 ) is so large as to substantially raise the surface 
temperature (the increase depending on various thermo- 
dynamical effects) during the impact periods. While this 
phenomenon has no bearing on the annihilation of our 
positrons, it can still lead to some (transitory) obser- 
vational effects, which may be investigated. Indeed, the 
depth at which positrons annihilate depends on the ther- 
modynamic conditions of the gas (density, temperature), 
and using the energy loss rate of positrons in a totally 
ionized medium (see Eq. 3 of Murphy et al. 2005 and ad- 
ditional references therein), we estimate that the depth 
at which positrons with initial kinetic energy of 660 keV 
annihilate is ~0.1 g/cm 2 for a temperature of 10 4 K and 
^0.06 g/cm 2 for a temperature of 10 6 K. Moreover, at 
this energy, the range of positron can vary also by a fac- 
tor of ~3.5 with the ionization fraction and the abundance 
(see Fig. 3 of Milne, The & Leising 1999). In any case, the 
range of 660 keV positron is lower than the mean free path 
of 511 keV photons (~6.6 g/cm 2 ). 

Even if a deformation of the surface does occur, the 
positrons still annihilate at a depth ~ A^ and the upward 
511 keV photons do escape the atmosphere. However, the 
impact of the jet produces a cavity which may reduce the 
escaping solid angle of the upward 511 keV photons. We 
take this effect into account in an approximate, geometric 
way. 

The jet pressure on the stellar atmosphere is Pj et — 
Lj e t/ (vj e tS) , where Lj et is the power of the jet, wj e t its 
velocity, and S is the impact surface on the secondary. As 
a first approximation, we consider that the jet stops in 
the secondary star's atmosphere when the jet pressure is 
equal to the atmospheric gas pressure. At the impact, the 
jet compresses the atmosphere, producing a cavity with 
a depth that depends on the pressure profile of the sec- 
ondary's atmosphere, the jet characteristics (power, veloc- 
ity and opening angle) and the binary separation. For typ- 
ical values of jet power Lj et ^10 36 erg s _1 , velocity Vj et ~ 
c/2 and opening angle 9j ct ~ 1°, we obtain Pj et ~ 10 5 
erg cm~ 3 , for a binary star separation of a = 2 R Q and a 
secondary mass of 0.5 M Q . Using the density and the tem- 
perature profile models of a 0.5 M secondary atmosphere 
(Guessoum & Jean 2002), we derived the atmospheric gas 
pressure as a function of the depth and estimated cavity 
depths (below the photosphere) of ~10 7 cm. 

The fraction of upward 511 keV photons that escape 
the cavity can be estimated by calculating the opening an- 



gle of the cavity (solid angle viewed from the center of the 
bottom of the cavity to the edge of the top of the cavity). 
Using the same jet-binary parameters (the radius of a 0.5 
M is ~0.7 Rq) and a jet impinging the surface perpen- 
dicularly, the radius of the jet impact on the secondary 
atmosphere is 7.9 x 10 s cm. With a depth of 10 7 cm, the 
opening angle is then ^178°, yielding a solid angle of 6.2 
sr, so the fraction of upward photons that escape the cav- 
ity is 98%. Considering the most pessimistic cases (small 
binary separation and jet impinging perpendicularly on 
the surface, i.e. high jet pressure), the fraction of upward 
511 keV photons escaping the cavity is always larger than 
70% (e.g. with a separation of a = 1.4 Q , the cavity depth 
is ^10 8 cm and the escaping fraction is ~77 %). Therefore, 
we neglect this jet-pressure effect on the secondary atmo- 
sphere in the modeling of the 511 keV flux presented in 
the next paragraphs. 

We also checked that the magnetic field of the sec- 
ondary does not affect the jet: we compared the jet pres- 
sure with the magnetic pressure Pb — B 2 /(4tt); taking a 
most extreme value for the stellar magnetic field B ^10 3 G 
we obtain Pg ~ 8 x 10 4 erg cm -3 which is lower than Pjet- 

The upward 511 keV photons that escape the atmo- 
sphere may enter the incoming jet. They traverse the jet 
if their mean free path is larger than the jet dimension. 
The minimum mean free path of 511 keV photons in the 
jet material is A > (neax) 1 , with n e the electron plus 
positron density in the jet and ot the Thompson cross- 
section. At the secondary surface, n e can be approximated 
by n e ~ 2N e + / (v ie tS). Assuming a close binary system 
(e.g. a = 2 R Q and R2 = 0.7 R©) to estimate the small- 
est impact area S and therefore the largest density n e , 
we obtain A > 6 R© (assuming Wj ct ~ c/2 and 6*j et ~ 
1°). This lower limit is already larger than the size of the 
XRB separation. Consequently, the jet itself can be con- 
sidered as quasi-transparent to those gamma photons and 
the general conclusion is that positrons annihilate near the 
surface, and about half the resulting 511 keV photons will 
be emitted upwards, towards the surface and the incoming 
jet. 

The binary system rotates, and the secondary star 
periodically crosses the jet, which slowly precesses. 
Consequently, the 511 keV photon radiation flux received 
by an observer in a given direction is a function of time, 
and its intensity and time-variability characteristics de- 
pend on the parameters of the XRB (period, inclination, 
distance, separation and radius of the secondary star) and 
of the jet (inclination with respect to the binary plane, 
opening angle and precession period). The flux at 511 keV 
can be calculated using the following expression: 

^siikevW = 2 x (1 - |/ Pa ) x N e+ x ^ (8) 

where fp s is the positronium fraction, N e + is the rate of 
ejection of positrons from the jet, and d is the distance of 
the binary system to the observer. The function I(t) is the 
fraction of 511 keV photons coming from the atmosphere 
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Fig. 2. Schematic view of the irradiation of the secondary 
star by a jet. The binary system frame is represented by 
the dashed axes. 



in the direction of the observer. This function takes into 
account the fraction of the jet's solid angle (which contains 
the positrons) that (i) intersects the secondary star surface 
and (ii) is visible by the observer. Its expression is given 
by: 

1 f 

I{t) = — / e *5iik e v<«>»ec*>.*> dil(<p, if)) (9) 

where A# is the atmospheric depth (in g cm -2 ) at which 
annihilation occurs, and AsukoV is the mean free path 
of 511 keV photons in the stellar atmosphere. Af2j ct is 
the solid angle of the jet; <p and if: gi ye the direction of 
the jet (see Fig. 0); and D{t) is the integration domain 
((p,tp £ D(t)), defined as the intersection of the irradiated 
area and the visible area in the direction of the observer 
(the set of directions ip and rp for which 6 < 90° providing 
the observer direction a and /?). The irradiation area (or 
hotspot) is delimited by the intersection of the compan- 
ion surface with the cone of angle 0j e t/2 around the jet 
direction cp$ et and ip^ ct (see Fig. [5J. In Equation 6 is 
also a function of a, f3, the XRB separation (a) and the 
secondary radius (i?2)- Since I(t) is expressed in the XRB 
frame that is rotating with respect to the observer, a, tpj e t 
and i/jjet vary with time, while (3, the angle between the 
binary system plane and the observer direction, is equal 
to 90° — i (i is the "official" inclination of the XRB). 

To illustrate the temporal variations of the annihila- 
tion flux, we calculated the 511 keV light curve for the 
case of a "prototypical misaligned microquasar" , one that 
has similar parameters (separation, inclination, primary 
and secondary masses. . . ) as GRO J1655-40, which is one 
of two XRB's for which independent information is avail- 
able on the jet orientation (Narayan & McClintock 2005). 
However, we arbitrarily assumed a jet inclination of ij Ct = 
9° with respect to the binary plane. With such an inclina- 
tion, the jets hit the companion star at about half of its 
upper hemisphere. We also assumed an opening angle #j e t 
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Fig. 3. Instantaneous 511 keV photon flux induced by the 
jet (red curve) and the counterjet (blue curve) impacts 
on the secondary of a virtual XRB (see text). The jet 
inclination with respect to the binary plane is taken to 
be ij e t = 9°, its opening angle is 0j ct = 1°; the precession 
period is 300 days (see text). The positron rate and energy 
in each jet is N e + — 10 41 s _1 and E — 660 keV (vj ot = 
0.9 c), respectively. 

= 1° and a precession period of 300 days which is the order 
of magnitude of the precession period of such an inclined 
jet (Kaufman Bernado, Romero & Mirabel 2002). 

Could there be some observational effects on the radio 
emission from the interaction of the jet with the surface 
of the star? In this regard, we first note that pairs in the 
jets will emit synchrotron radiation for most (say 90 %) 
of their transit (ejection) time, while the rest of the time 
they spend depositing their energies in the companion's 
atmosphere; this energy is lost first and foremost in the 
form of heat and then, at the end, as gamma rays, which 
get downgraded as X-rays and again as heat. It is thus not 
expected that much radio emission will come out of the 
jet's collision with the companion star's atmosphere, but 
this issue could be investigated further in the future. 

As for GRO J1655-40, this "prototypical" microquasar 
was assumed to be at a distance of 3.2 kpc from Earth, 
to have a period of 2.62 days and to be composed of a 
compact object with a mass of 6.8 M Q separated by 16.6 
R© from a companion star with a radius of 5 R0. The 
inclination of the binary system is i = 70°. Figure 01 shows 
the instantaneous 511 keV flux emitted by such a micro- 
quasar as a function of time. Details of this light curve are 
presented in Figure0]for two different dates: (a) when the 
direction of the jets lies approximately along the line of 
sight, and in (b) when the direction of the jets is approxi- 
mately perpendicular to the line of sight. The geometrical 
configurations of the jets and the binary system in these 
two cases are shown in Fig. |SJ 

Figure [S] shows the time series of the 511 keV flux 
averaged per period of the XRB. The flux is minimum 
when the jet direction is perpendicular to the observer 
direction, as shown in Fig. 0d. In this configuration, the 
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Fig. 4. Details of the instantaneous 511 keV photon flux 
induced by the jet (red curve) and the counterjet (blue 
curve) impacts on the secondary of a virtual XRB (see 
caption of Fig.[3Jl. (a) The jet direction is parallel to the 
observer direction and only the hotspot induced by the 
first jet is visible; (b) the jet direction is perpendicular to 
the observer direction and both hotspots induced by the 
jet and the counterjet are visible. The geometrical config- 
urations of the jets and the binary system in these two 
cases are shown in Fig. [S] 



inclination 9 of the direction of the photons emitted by the 
fraction of the hotspot that is visible can be close to 90° 
during the crossing of the secondary. Consequently, a large 
fraction of these photons are absorbed in the atmosphere 
of the secondary (see Eq. [5] and Fig. HJd). Note that the 
average flux is ~10 times lower than the maximum of the 
instantaneous flux since the companion star crosses the 
jets during a fraction ~ 10% of the period. 

The spectral shape of the annihilation radiation emit- 
ted by this virtual microquasar was calculated accordingly 
to the model of Guessoum, Jean & Gillard (2005), which 
provides annihilation spectra as a function of the physical 
conditions of the gas in which the annihilation occurs. We 
assumed that positrons annihilate in a half-ionized and 
half- neutral plasma with a temperature of 8000 K. The 
rotation of the binary system leads to a shift in the cen- 
troid of the 511 keV line. This Dopplcr shift is a function 
of time and depends on the observer direction. We cal- 
culated this spectral shift using the method described in 
Jean & Guessoum (2001). Figure shows the spectra: in- 
tegrated over the precession period (case 1); integrated 
over one XRB period at t = 100 days (case 2, see Fig.[3J); 
integrated over one XRB period at t = 200 days (case 3). 
The line is blue-shifted by 0.3 keV in case 2 because most 
of the annihilation emission comes from the jet impact 
on the secondary which is moving toward the observer. 
In case 3, most of the annihilation emission comes from 
the jet impact on the secondary which is receding from 
the observer, thus producing a line with a red-shift of 0.3 
keV. 




Fig. 5. Schematic views of the jets (blue) arising from the 
compact object (black), the secondary star (red) presented 
at four positions in its orbit, and the annihilation hotspot 
(yellow) emitted from the secondary star atmosphere, pre- 
sented at the dates ~ 2 days (a) and ~ 75 days (b). The 
directions of motion of the secondary star and jets axis are 
represented by arrows. The two configurations correspond 
approximately to the cases a and b of Fig. 21 
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Fig. 6. Average 511 keV photon flux per XRB period, 
induced by the jet impacts on the secondary of an XRB. 
The parameters of the jets are identical to those of Fig.|3J 



The results presented in this section show that the 
measurements of the temporal and spectral profiles of the 
annihilation radiation emitted by microquasars with mis- 
aligned jets are useful to probe the characteristics of the 
jet. The maximum 511 keV instantaneous flux provides 
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Fig. 7. Spectra of the annihilation radiation emitted by 
GRO J1655-40. The parameters of the jets are identical 
to those of Fig. [3] We assumed that the 50% of positrons 
annihilates with the ionized component of the atmosphere 
and 50% with the neutral component, case 1: total spec- 
trum integrated over the precession period; case 2: spec- 
trum integrated over one period at the date 100 days; 
case 3: spectrum integrated over one period at the date 
200 days. 

the rate of positrons in the jet (e.g. Fig- EJ) - The duration 
of the 511 keV emission in one XRB period allows one to 
estimate the opening angle of the jet (e.g. Fig. while 
the long-term variation of the average flux allows for a de- 
termination of the precession period of the jet (e.g. Figj|jJ). 
The spectral shape of the 511 keV line, do, in principle, 
allow one to infer the physical conditions of the plasma in 
which positrons annihilate, as Jean et al. (2006) demon- 
strated for the SPI data of positron annihilation in the 
ISM, although one must realize that in more complex sit- 
uations a host of conditions and effects can render the 
extraction of physical information very complicated if not 
impossible. 

Finally, in Table [21 we present our predicted maximum 
fluxes for the confirmed misaligned microquasars as well 
as for potential candidates, i.e. microquasars that may 
be at least occasionally misaligned and are not too far 
away for detection by future instruments. We must note 
that the values given in this table are for "average jets" ; 
"strong-flare jets" could be 10 times more powerful, thus 
producing 10 times more positrons and thus resulting in 
an annihilation flux 10 times higher (than what is given in 
the table); on the other hand, "weak episode jets" would 
give positron rates and annihilation fluxes about 10 times 
lower. Indeed, the annihilation flux is directly proportional 
to the jet power. 

6. Summary 

The origin of the Galactic positron annihilation radiation 
is a major issue in high energy astrophysics at present. In 
this work we assess the potential of Galactic microquasars 



(XRBs that exhibit jets in an intermittent way) as sources 
of positrons. 

Among the rapidly growing body of data on micro- 
quasars, we present in Sec. 2 the main features, those of 
relevance to our study. The correlation between the power 
of the jet and the X-ray luminosity of the compact ob- 
ject is the most important of these features. We stress, in 
particular, that the ratio between the two, often assumed 
to be small, is highly uncertain; moreover, the content of 
the jets, leptonic or baryonic (i.e. electron-positron pairs 
vs. protons and pions), is unknown at present. Some of 
the implications of a baryonic content have been explored 
elsewhere (Butt, Maccarone & Prantzos 2003, Romero & 
Orelana 2005) ; we explore here the consequences of a lep- 
tonic content. 

In Sec. 3 we evaluate the rate of positron ejection by 
the microquasar jets, based on various models proposed 
in the literature, but also on simple energetic arguments 
relating the total power of the positrons in the jets to the 
estimated total X-ray luminosity of the "low luminosity - 
hard spectrum" galactic LMXRBs. We find that a value of 
10 41 e + s _1 could be considered as a "canonical average", 
albeit with large (and difficult to evaluate) uncertainties. 

In Sec. 4 we estimate the total positron production 
rate by the collective emission of Galactic microquasars 
and find it smaller (about one third) than what is in- 
ferred from SPI /INTEGRAL observations of 511 keV ra- 
diation. The spatial morphology of the corresponding flux 
received on Earth depends on the assumed large scale dis- 
tribution of Galactic microquasars, which is very poorly 
known at present; we find that the distribution of the cur- 
rently available (incomplete) sample appears encouraging, 
in that respect. (We should also note that a large enough 
B/D ratio would be obtained with any source distributions 
that congregate in the bulge and/or suppress the annihi- 
lation of their positrons in the disk.) Finally, we constrain 
the rate of production of positrons by microquasars on the 
basis of the SPI flux results: we find that the SPI limits 
would be violated if positrons are produced at steady rates 
greater than ~3x 10 41 e+s^ 1 . 

Finally, in Sec. 5 we present an investigation of the 
physics of "misaligned" microquasars, where the jet im- 
pinges on the companion star. We find that about half 
of the annihilation photons should, in general, escape the 
stellar atmosphere. We compute fluxes for a list of sources 
(see Table |2Jl, a few of them known to be misaligned, and 
some of the others may turn out to be so as well. The 
expected annihilation fluxes are found to be low, how- 
ever, and difficult to measure with the spectrometer SPI 
due to the distance of most of the "misaligned" micro- 
quasars. Future instruments such as a gamma-ray lens 
(von Ballmoos et al. 2004) or a Compton telescope (Boggs 
& Jean 2001) would have the sensitivity required to mea- 
sure these fluxes. Due to the rotation of the binary system, 
the secondary star periodically crosses the jet, leading to 
a "periodic" and a Doppler-shifted annihilation emission. 
Due to the precession of the jet, we expect to observe 
a long-term modulation of the emission. A detection of a 
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Table 2. Predicted maximum 511 keV emission fluxes for confirmed or potentially misaligned nearby microquasars. 
Fluxes were calculated for "Average Jet", that is using a positron rate of 10 41 e + s _1 . 



Source 


Jet 


Distance 


Flux for "Average Jet" 






(kpc) 


(ktWs- 1 ) 


GRO J 1655-40 


Misaligned 


3.2 


0.49 


V4641 Sgr 


Misaligned 


9.6 


0.05 


XTE J1550-564 


Misaligned? 


5.3 


0.18 


X1E J1118+480 




1.8 


1.53 


LS I +61° 303 




2 


1.24 


Cyg X-l 




2.5 


0.79 


Sco X-l 




2.8 


0.63 


LS 5039 




2.9 


0.59 


SS 433 




4.8 


0.22 


GRS 1758-258 




8.5 


0.07 


Cyg X-3 




9 


0.06 


Cir X-I2 




10 


0.05 


GRS 1915+105 2 




12.5 


0.03 


GX 339-4 




> 8 


< 0.08 


XTE J1748-288 




> 8 


< 0.08 


IE 1740.7-2942 




8.5 


0.07 


IGR J17091-3624 




8.5? 


0.07 


IGR J17303-0601 




8.5? 


0.07 


IGR J17464-3213 




8.5? 


0.07 


IGR J 18406-0539 




8.5? 


0.07 



point source at 511 keV from a "misaligned" microquasar, 
with such temporal and spectral signatures would con- 
firm the presence of positrons in the jet. Moreover, such a 
measurement would allow for a determination of the rate 
of positrons channeled by the jet and the characteristics 
(opening angle, precession period) of the jet. 

In summary, we have shown in this study that micro- 
quasar jets may constitute important sources of annihila- 
tion radiation, both diffuse (for most of them) or point 
sources (in the case of misaligned ones). 
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